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ABSTRACT 


.-'Thin  study  covers  the  development  of  a  set  of  computer  programs  to 
perform  flutter  analysis  by  the  collocation  method.  While  this  method 
has  been  known  for  some  time,  only  recently  have  advances  in  computer 
technology  made  the  method  technically  and  economically  feasible.  The 
ingredients  of  a  collocation  flutter  analysis  are  (l)  a  flexibility 
matrix, 12)  aerodynamic  influence  coefficient  matrix,  and'3)  an  eigenvalue 
solution.  This  study  is  presented  in  four  volumes.  N  Volume  I  contains  a 
general  program  discussion./  -Volume  II  contains  t)*e  program  FLUESC  which 
calculates  the  flexibility  matrix.  Volume  III  contains  a  set  of  three 
programs  to  calculate  aerodynamic  influence  coefficients  for  subsonic, 
transonic,  and  supersonic  flight  regimes.  Volume  IV  contains  the  program 
COFA  which  sets  up  and  solves  the  flutter  eigenvalue  matrix. 
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FOREWORD 


This  report  covers  the  research  conducted  by  the  Missile  Systems 
Division  of  Hughes  Aircraft  Company,  Canoga  Park,  California,  under 
Contract  Number  00019-68-C-0247. 

This  work  was  performed  to  develop  a  comprehensive  set  of  computer 
programs  that  perform  flutter  analysis  using  the  collocation  method.  A 
set  of  three  computer  programs  to  calculate  unsteady  aerodynamic  influence 
coefficients  for  subsonic,  transonic,  and  supersonic  lifting  surfaces 
with  interaction  between  the  wing  and  a  downstream  .control  surface  lying 
in  the  wake  of  the  wing  are  presented.  A  program  to  calculate  structural 
influence  coefficients  is  presented,  and  a  program  to  set  up  and  solve 
the  flutter  eigen  matrix  is  presented. 

The  Program  Manager  for  Hughes  Aircraft  Company  was  R.  J.  Oedy. 

Dr.  V.  Weingarten  assisted  in  the  development  of  the  structural  analysis 
program,  and  Dr.  W,  P.  Rodden  assisted  in  the  development  of  the  aero¬ 
dynamic  programs  and  flutter  analysis  programs. 
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REFERENCES 


PART  I  -  INTRODUCTION  | 

a 


Collocation  is  a  tehn  used  to  describe  a  mathematical  procedure  for 
(among  other  things)  solving  integral  equations.  The  word  literally  means 
to  co-locate,  i.e.,  an  approximate  solution  is  co-located  to  satisfy  the 
equation  at  a  specified  discrete  number  of  locations  in  the  space  (the 
collocation  points).  When  the  integral  equation  for  vibration  analysis  is 
solved  by  the  collocation  method  it  is  generally  referred  to  as  the  lumped- 
parameter  method  since  the  properties  of  the  continuum  are  assumed  to  be 
concentrated  ("lumped")  at  the  collocation  points.  The  collocation  method 
leads  to  a  formulation  of  the  flutter  problem  similar  to  that  of  the 
vibration  problem  where,  in  addition  to  the  structural  influence  coefficients 
(SIC’s)  and  mass  matrix,  aerodynamic  influence  coefficients  (AIC's)  are 
required. 

The  vibration  problem  can  be  large  if  many  degrees  of  freedom  are 
necessary  to  model  the  vibrating  system,  and  the  calculation  of  the  modes 
and  frequencies  requires  solving  a  large  eigenvalue  problem.  However,  since 
both  the  SIC’s  and  mass  matrix  are  real  symmetric  matrices,  special  techniques 
can  be  utilized  to  find  the  modes  and  frequencies  rapidly  (and,  hence, 
economically).  The  corresponding  formulation  of  the  flutter  problem  does  not 
possess  any  special  properties  that  permit  convenient  determination  of  the 
modes,  frequencies,  and  flutter  velocities.  For  this  reason  a  series  or 
modal  solution  has  traditionally  been  employed  to  reduce  the  size  of  the 
flutter  problem  and  the  corresponding  computing  time. 

Present  day  computer  capacities  nnd  speeds,  and  advances  in  eigenvalue 
methods  however,  no  longer  limit  flutter  analysis  to  modal  methods  and  the 
mi  vantages  ol'  a  direct  solution  can  be  realized.  One  advantage  is  that 
question:;  of  convergence  of  the  modal  solution  are  avoided.  Another  advantage 
is  thnl  the  flexibility,  mnar.  nnd  aerodynamic  data  are  kept  separate  so  that 
n  change  in  nny  one  set  cnri  be  made  independently  of  the  other  two.  Finally, 
the  flutter  nnnlyr.is  can  ho  done  without  performing  a  vibration  analysis.  This 
is  n  particular  advantage  when  variations  in  mass  or  stiffness  distribution  must 
be  Investigated,  lixamples  include  control  surface  mass  balance  optimization 
studies  or  weight  reduction  due  to  burning  fuel. 
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PART  II 


GENERAL  DISCUSSION  OF  FLUTTER 


A.  The  Mechanism  of  Flutter 

Flutter  Is  a  self-excited  vibration  of  a  flexible  body  in  a  fluid  stream. 
The  flutter  speed  is  defined  as  the  minimum  speed  at  vhich  the  aeroelastic 
system  will  sustain  a  neutrally  stable  oscillation.  At  lover  speeds  the 
motion  will  be  damped,  and  at  higher  speeds,  or  at  least  in  a  range  of  higher 
speeds,  the  motion  will  be  divergent.  Usually  only  a  small  increase  beyond 
the  flutter  speed  is  required  to  produce  a  rapid  divergence  of  such  violence 
that  a  complete  structural  failure  results  in  only  a  few  cycles  of  motion. 

At  the  flutter  speed,  a  critcial  phasing  between  the  motion  and  the 
loading  permits  extraction  of  an  amount  of  energy  from  the  air  stream  equal 
to  that  dissipated  by  internal  damping  during  each  cycle  and  thereby  sustains 
a  neutrally  stabla  periodic  motion.  The  simplest  type  of  flutter  occurs 
when  a  simple  motion  induces  an  aerodynamic  force  having  a  component  in  the 
direction  of  the  motion  and  in  phase  with  the  velocity.  This  is  described 
as  a  case  of  negative  aerodynamic  damping  or  single  degree  of  freedom 
flutter.  The  term  classical  flutter  is  used  to  denote  the  more  complicated 
instability  that  typically  arises  from  a  critical  coupling  of  two  or  more 
modes  of  motion,  each  of  which  is  usually  stable  by  itself.  It  is  the 
classical  flutter  mechanism  with  which  we  are  primarily  concerned  because 
its  generality  includes  the  single  degree  of  freedom  instability  as  a 
special  case. 

B.  The  Ingredients  for  Flutter  Analysis 

Since  flutter  is  a  dynamic  aerolastic  instability  the  forces  involved 
arc  structural,  aerodynamic,  and  inertial.  In  the  present  lumped-parameter 
or  collocation  formulation  these  appear  as  matrices  of  structural  influence 
coefficients  (.'SIC's),  aerodynamic  influence  coefficients  (AIC’s),  and  masses. 
The  combination  of  the  various  matrices  to  formulate  the  flutter  stability 
problem  lends  to  a  non-Hermitian  eigenvalue  problem.  Since  this  is  the  most 
general  type  of  eigenvalue  problem,  methods  for  its  solution  are  quite  limited. 
The  method  that  has  been  employed  is  the  power  method,  or  as  it  is  sometimes 
called,  the  matrix  iteration  method. 
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PART  III 


DISCUSSION  OF  COFA  PROGRAM 


A.  The  Flutter  Analysis 

1 

I.  The  Technical  Approach  - 

The  collocation  method  of  flutter  analysis  had  been  known  for  many 
years  but  it  was  not  demonstrated  to  be  computationally  feasible  until  1956 
in  Reference  1.  The  method  was  extended  to  include  an  arbitrary  number  of 
rigid  body  degrees  of  freedom  in  Reference  2.  A  complete  discussion  of  the 
general  flutter  analysis  is  given  in  the  Flutter  Analysis  Report.  It  is 
sufficient  for  the  introductory  purposes  of  this  report  to  consider  the 
formulation  of  the  flutter  problem  for  a  restrained  vehicle. 

The  basic  equation  is  the  relationship  between  the  deflections  f  h] 
to  the  forces  through  the  structural  influence  coefficients  (SIC’ a) 

r  • 

,  a  t  . 

i  J 


fh)  -  M  W 


(2.1) 


The  force  matrix  must  include  both  the  inertial  end  aerodynamic 


forces. 


F  » 


r  -i  C"l  ^  1 

lM]  W  *  iFai 


(3.2) 


where  M  is  the  mass  matrix  corresponding  to  the  chosen  set  of  control 

w  * 

points.  We  define  a  complex  matrix  of  oscillatory  aerodynamic  influence 
coefficients  (AIC's),  j  C.1, by  the  equation 


jF„)  ’  *r  8  [°h]  (h] 


(3.3) 


where  p  ir.  the  density,  us  is  the  frequency,  br  is  the  reference  semi-chord, 
and  s  Jr.  the  reference  serai-span.  Combining  Eqs.  (3.2)and  (i.3)  for  harmonic 
motion,  we  find  the  total  force  to  be 


M  -  “2  (K  *  pbr2”  (v)  H 


(3.4) 
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If  we  substitute  Eo.  (3,4)into  Eq.  (3,1) and  divide  the  static  SIC's 
by  (l  +  Ik)  where  g  is  the  artificial  structural  damping  necessary  to 
sustain  the  assumed  harmonic  motion,  we  obtain  the  matrix  aquation  for 
f lutter 

Ww  [*](["]  K])  W 

2q.(3.5)  can  be  solved  for  the  complex  mode  shape,  £h]  ,  and  the  complex 

eigenvalue,  (1  +  ig)/<v^  by  complex  matrix  iteration  (the  power  method). 
From  the  eigenvalue 


XR  +  ikj 

(3.6a) 

-  (1  +  ig)/«2 

(3.6b) 

we  obtain  the  flutter  frequency 

w  *=  1/  ^  X  R 

(3,7) 

the  required  artificial  damping 

6  ~ 

(3.8) 

and,  since  i.he  AIC's  required  the  assumption  of  a  reference  reduced  frequency, 
kr  for  its  calculation,  we  find  the  velocity 

U  =ubykr  (3,9) 

From  a  series  of  solutions  for  different  reduced  frequencies,  a 
required  damping  velocity  curve  can  be  constructed  for  a  specific  altitude. 

The  flutter  velocity  corresponds  to  the  lowest  velocity  at  which  the 
required  damping  is  equal  to  the  actual  structural  damping  of  the  system. 
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Program  Capabilities 

The  COFA  program  has  the  limitations  of  analyzing  a  system  with  a  maximum 
of  49  collocation  stations.  In  the  use  of  these  49  stations,  there  is  a  great 
deal  of  flexibility.  Entire  vehicles  may  be  analyzed  as  single  entities  or  as 
conglomerates  of  components  with  single  or  tanden  aerodynamic  surfaces.  Both 
may  be  analyzed.  When  the  components  are  used  to  synthesize  a  structure,  the 
attachment  between  the  components  must  be  statically  determinant.  The 
components  may  be  either  flexible  or  rigid.  The  coupling  of  flexible  components 
is  demonstrated  in  the  example  problem  of  this  volume.  The  coupling  of  flexible 
to  rigid  components  is  demonstrated  in  the  example  problems  of  Volume  IV.  The 
systems  analyzed  may  be  cantilevered  or  free  in  space.  Up  to  six  rigid  body 
modes  may  be  used  in  the  analysis  of  systems  free  in  space.  In  addition,  there 
is  the  capability  of  using  a  free- free  boundary  condition  built  into  the  program 
for  freeing  cantilevered  or  fixed  systems. 

There  are  two  classes  of  tandem  surface  problems;  they  are  aerodynamically 
non- interfering  surfaces  and  aerodynamically  interfering  surfaces.  For  non¬ 
interfering  surfaces  the  structural  system  may  be  analyzed  as  a  single  entity 
or  a  conglomerate  of  components.  The  most  usual  way  is  to  use  the  component 
synthesis  technique.  When  the  component  technique  is  used,  the  AIC,  SIC,  and 
mass  matrices  are  partitioned  with  a  partition  corresponding  to  each  structural 
component.  This  results  in  dealing  with  smaller  order  matrices,  and  the 
analysis  is  straightforward  as  presented  in  Volume  IV.  The  only  precaution 
is  that  the  mass,  flexibility,  and  aerodynamic  matrix  partitions  be  entered 
serially  into  the  computer  program  in  an  identical  manner,  so  that  the 
iniormation  is  properly  located  in  the  overall  dynamical  matrix,  i.e.  fuselage 
first;  forward  surface  second;  and  aft  surface  third;  etc.  For  interfering 
mir I  aces,  the  analysis  must  be  conducted  as  a  single  structural  entity;  for 
Inter  Luring  surfaces  the  following  procedure  must  be  followed: 
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If  fuselage  aerodynamics  are  obtained  from  a  different  source,  they 
may  be  properly  entered  into  the  analysis  in  the  upper  left  partition.  The 
off  diagonal  partitions  are  for  wing-fuaelage  interference  effects.  The 
details  of  using  the  COFA  Program  may  be  found  in  Volume  IV. 


B.  The  Calculation  of  Structural  Influence  Coefficients 
1.  The  Technical  Approach  - 

There  are  many  methods  by  which  structural  influence  coefficients 
may  be  calculated.  The  SIC's  calculated  for  this  program  use  the  direct 
stiffness  method.  The  direct  stiffness  method  was  chosen  as  it  is  readily 
adapted  to  highly  redundant  structures  and  is  especially  suited  to  the  use 
of  high  speed  digital  computers.  The  structure  is  regarded  as  an  assemblage 
of  beams  and  plates.  The  stiffness  matrix  for  the  entire  structure  is 
computed  by  the  simple  summation  of  the  stiffness  matrices  of  the  elements 
of  the  structure.  Finally,  the  matrix  of  SIC's  is  obtained  by  the  inversion 
of  the  stiffness  matrix. 


For  *  given  idealized  structure,  the  analysis  of  stresses  and 
deflections  due  to  a  given  system  of  loads  is  a  purely  mathematical 
problem.  Two  conditions  must  be  satisfied  in  the  analysis:  1)  the 
forces  developed  in  the  members  must  be  in  equilibrium  and,  2)  the  deform¬ 
ation  of  the  members  must  be  compatible  -  l.e.,  consistent  vith  each  other 
and  the  boundary  conditions.  In  addition,  the  forces  and  deflections  in 
each  member  must  be  related  in  accordance  with  the  stress-strain  relationship 
assumed  for  the  material.  These  conditions  are  satisfied  by  the  direct 
stiffness  method  in  the  following  way.  The  displacements  of  the  joints  of 
the  structure  are  considered  as  unknown  quantities.  Thus,  an  infinite 
number  of  systems  of  mutually  compatible  deformations  in  the  members  are 
possible;  the  correct  pattern  of  displacements  is  the  one  for  which  the 
equations  of  equilibrium  are  satisfied.  The  concept  of  static  detenainate- 
ness  or  indeterminateness  is  irrelevant  when  the  analysis  is  considered 
from  this  viewpoint.  This  approach  is  the  basis  for  many  relaxation  type 
analysis  (such  as  moment  distribution)  and  is  used  extensively  in  the 
analysis  of  complex  aircraft  structures. 

2.  Program  Capabilities  - 

The  program  FLUENC  calculates  structural  influence  coefficients  for 
systems  with  up  to  50  mass  points.  The  mass  points  may  be  connected  by  any 
combination  of  beams  or  plates.  Each  member  may  have  its  own  material 
property.  The  program  also  calculates  the  natural  frequencies  and  mode 
shapes  of  the  structure.  The  program  can  calculate  the  modes  and  frequencies 
for  structures  fixed  or  free  in  space.  Caution  should  be  exercised  not  to 
use  the  flexibility  matrix  of  structures  free  in  space  as  the  flexibility 
matrix  of  such  a  structure  cannot  be  directly  calculated.  The  flexibility 
matrix  of  a  free  structure  can  be  obtained  by  l)  arbitrarily  fixing  the 
structure,  2)  calculating  the  flexibility  matrix  for  the  fixed  structure, 

3)  enter  the  flexibility  matrix  of  the  fixed  structure  into  the  COFA  program 
and  use  the  free-free  boundary  condition  to  obtain  the  flexibility  matrix 
of  the  free-free  atructure.  The  details  of  using  the  program  FLUENC  may 
bo  found  in  Volume  II. 


C*  The  Unsteady  Aerodynamic  Lifting  Surface  Theories 

1.  The  Subsonic  Kernel  Function  Procedure 
3)  The  Technical  Approach  - 

The  linearized  equation  for  the  perturbation  potential  <j>  for 
a  body  immersed  in  a  compressible,  imriscid,  perfect  gas,  when  we  assume 
tr.e  flow  to  be  isentropic  and  irrotational,  is 


(i  -  m2)  y  4> 


2  2 

3  4>  0  4> 

+  -  +  - 

3  y2  a  ,2 


+  tiL 


2 

a  4> 
ax  at 


(3.10) 


in  which  M  is  the  Mach  number,  aQ  is  the  ambient  speed  of  sound,  and  x,y.,  z 
are  the  Cartesian  coordinates  of  the  flow  field.  The  linearized  perturbation 
pressure  coefficient  is  related  to  the  potential  through 


S  ■ 


where  V  is  the  free-stream  velocity.  The  problem  in  lifting  surface  theory 
is  to  find  the  difference  in  the  pressure  coefficients  between  the  upper 
and  lower  surfaces,  i.e.,  the  lifting  pressure  coefficient.  The  evaluation 
of  Eo. (311)  requires  the  solution  of  Eq.  (3. 10). The  differential  equation  (3.10) 
must  be  solved  subject  to  certain  boundary  conditions.  There  can  be  no  flow 
through  the  surface,  i.e.,  at  the  surface  the  normal  velocity  of  the  fluid 
must  ocmal  the  normal  velocity  of  the  surface  (the  downwash  velocity). 
Additional  conditions  at  subsonic  speeds  require  that  the  pressure  must 
remain  finite  at  the  trailing  edge  (the  Kutta  condition)  and  the  potential 
must  vanish  at  an  infinite  distance  from  the  surface.  The  additional 
conditions  st  supersonic  speeds  are  that  no  disturbance  can  propagate  upstream 
of  the  lending  edge,  and  thnt  the  disturbances  must  propagate  from  the  surface 
nc  outgoing  waves.  At  transonic  speeds  a  combination  of  the  subsonic  and 
supersonic  boundary  conditions  are  required. 


(3.11) 
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The  foregoing  general  review  of  lifting  surface  theory  leads 
us  to  consider  specific  Mach  number  regimes.  The  first  numerical  method 
for  solving  the  subsonic  pressure  distribution  for  planar  lifting  surfaces 
undergoing  simple  harmonic  motion  was  developed  by  Watkins,  Woolston,  and 
Cunningham  (Reference  3)*  The  approach  taken  was  to  recast  the  differ¬ 
ential  equation  for  the  velocity  potential  into  an  integral  equation 
relating  the  surface  downwash  to  the  lifting  pressure  coefficient  (the 
acceleration  potential).  The  solution  was  carried  out  by  assuming  a  series 
form  for  the  pressure  coefficient  based  on  known  exact  solutions,  e.g. , 
the  two-dimensional  solution  for  chordwise  pressure,  and  the  elliptical 
spanwise  loading  on  an  elliptical  wing.  From  this  series  of  pressure  modes 
the  problem  reduced  to  solving  for  the  amplitude  of  each  pressure  mode. 

A  collocation  method  was  used  in  which  the  downwash  was  matched  at  a  discrete 
number  of  points.  The  numerical  difficulties  are  numerous  since  the  integral 
equation  is  singular  and  the  many  integrals  of  the  pressure  modes  require 
extreme  accuracy.  Since  the  publication  of  Reference  3,  many  investigators 
have  refined  the  procedures,  notably  Hsu  (Reference  4),  Vivian  and  Andrew 
(Reference  p)>  and  Rove  (Reference  6).  The  present  study  is  an  extension 
of  References  7  and  8  to  obtain  AIC's,  thereby  eliminating  the  limitations 
introduced  by  the  requirements  for  polynomial  surface  fitting  in  References 
7  and  8.  The  details  of  the  AIC's  and  the  numerical  problems  of  the  sub¬ 
sonic  lifting  surface  theory  are  discussed  in  the  Unsteady  Aerodynamics 
Report,  Part  III  and  Part  IV,  Section  A,  respectively. 
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b)  Program  Capabilities  - 

The  Subsonic  AIC  program  calculates  AIC's  for  systems  with  up 
to  40  collocation  stations.  The  systems  analyzed  may  be  a  single  surface 
or  tandem  surfaces.  Up  to  20  collocation  stations  per  surface  may  be  used. 

When  analyzing  single  surfaces,  a  dummy  trailing  surface  must  be  input  and 
isolated  from  the  leading  surface.  The  following  single  surface  planforms 
may  be  analyzed:  rectangle,  trapezoidal,  delta,  cropped  delta  and  cropped 
trapezoid.  The  trailing  edge  of  the  surface  must  be  straight  and  normal  to 
the  direction  aerodynamic  flow.  When  analyzing  tandem  surfaces,  the  forward 
surface  may  be  anyone  of  the  single  surfaces  previously  defined;  the  aft 
surface  must  be  a  rectangular  surface  with  the  spanwise  dimension  equal  to 
the  spanwise  dimension  of  the  wing. 

The  program  automatically  sets  the  number  of  chordwise  and  spanwise 
pressure  series  terms  equal  to  the  number  of  chordwise  and  spanwise  collocation 
station.  This  maximizes  the  number  of  pressure  series  terms  used  and 
increases  computation  accuracy.  For  best  accuracy  a  minimum  of  three  chord- 
wise  collocation  stations  should  be  used.  The  number  of  kernel  function 
chordwise  integration  stations  is  a  variable  in  the  program,  and  is  a 
function  of  the  number  of  chordwise  collocation  stations.  (The  number  of 
integration  stations  equals  factor  times  the  number  of  chordwise  collocation 
stations.)  The  factor,  NI0NCX,  which  is  input  data,  should  be  such  that  the 
minimum  number  of  integration  stations  is  15,  The  maximum  number  of 
integration  stations  per  surface  is  40. 

When  analyzing  single  surfaces,  the  dummy  surface  is  required 
to  have  two  chordwise  collocation  stations.  The  number  of  spanwise  collo¬ 
cation  stations  on  the  dummy  surface  is  established  internally  in  the  program 
and  is  set  equal  to  the  number  of  spanwise  collocation  stations  used  in  the 
analysis  of  the  single  surface.  The  details  of  using  the  Subsonic  AIC 
program  may  be  found  in  Volume  III, 


11 


2.  The  Transonic  Box  Method 


a)  The  Technical  Approach  - 

The  basic  development  of  the  transonic  box  method  vas  made  by 
Roderoich  and  Andrew  (Reference  9)»  References  7  and  8  include  extensions 
of  Reference  9  to  interfering  wing-tail  configurations.  The  basic  formu¬ 
lation  of  the  transonic  potential  equation,  pressure  coefficient,  and 
boundary  conditions  has  been  included  in  the  introductory  discussion  of 
the  technical  approach  to  the  subsonic  problem  above  in  Section  A,  and  the 
details  of  the  AIC’s  and  the  numerical  problems  of  the  transonic  box 
method  are  discussed  in  the  Unsteady  Aerodynamics  Report,  Part  III,  and 
Part  V,  Section  A,  respectively.  For  the  purposes  of  this  introductory 
summary  it  is  sufficient  to  outline  the  basic  procedure. 

The  transonic  box  method  is  designed  to  calculate  the  unsteady 
potentials  from  which  the  pressure  distributions  may  be  obtained  for 
arbitrary  modes  of  surface  motion.  The  method  used  was  suggested  by  the 
successes  of  the  supersonic  box  methods  (References  10  and  11).  The 
potential  is  generated  by  a  doublet  distribution  rather  than  by  a  source 
distribution  because  the  latter  method  would  involve  diaphragm  regions  of 
infinite  extent,  whereas  the  doublet  distribution  is  confined  to  the  wing 
and  its  wake.  As  with  the  subsonic  problem,  the  differential  equation 
is  rewritten  as  an  integral  equation.  The  integral  equation  is  approximated 
numerically  by  a  matrix  equation  so  that  the  basic  step  in  the  box  method 
is  the  solution  of  the  system  of  simultaneous  equations  which  determine 
a  set  of  values  of  potential  on  the  surface  from  a  corresponding  array  of 
uprwash  values.  The  solution  procedure  obtains  the  velocity  potential  over 
the  surface  one  spanwise  row  of  boxes  at  a  time  urttil  the  trailing  edge  row 
is  completed.  The  numerical  complexity  is  not  increased,  however,  by  a 
large  number  of  box  rows  over  the  configuration  because  the  influence 
coming  from  more  than  15  rows  away  has  been  found  to  be  negligible.  The 
results  are  valid  for  high  reduced  frequency  k,  such  that  k  »  J  M  -1  1 
where  M  is  the  Mach  number. 
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b)  Program  Capabilities  - 

The  Transonic  AIC  program  calculates  AIC’s  for  systems  with  up  to 
40  collocation  stations.  The  systems  analyzed  may  be  a  single  surface  or 
tandem  surfaces.  Up  to  20  collocation  stations  per  surface  may  be  used. 

The  following  single  surface  planforms  may  be  analyzed:  rectangular, 
trapezoidal,  delta,  cropped  delta,  and  cropped  trapezoid.  It  is  not  necessary 
to  enter  a  second  surface  into  the  program  when  analyzing  a  single  surface. 

The  trailing  edge  of  the  surface  must  be  straight  and  normal  to  the  aero¬ 
dynamic  flow  direction.  When  analyzing  tandem  surfaces,  the  forward  surface 
may  be  any  one  of  the  single  surfaces  previously  defined;  the  aft  surface 
must  be  a  rectangular  surface  with  the  spanwise  dimension  equal  to  the 
spanwise  dimension  of  the  wing. 


The  maximum  allowable  number  of  Mach  boxes  is  45.  The  layout 
of  the  boxes  on  the  surface  is  determined  in  the  program,  and  is  a  function 
of  the  number  of  chordwise  boxes.  The  Mach  boxes  are  square,  and  can  be 
located  on  the  surface  by  placing  the  center  of  one  box  at  the  intersection 
of  the  root  chord  and  the  trailing  edfce  of  the  forward  surface.  The  length 
of  the  side  of  a  Mach  box  is  approximately  equal  to  the  root  chord  length 
divided  by  the  number  of  chordwise  boxes  chooaen.  It  is  not  necessary  to 
layout  the  Mach  boxes,  as  this  is  done  by  the  program,  but  this  knowledge 
is  helpful  in  choosing  the  optimum  number  of  chordwise  Mach  boxes.  In 
general,  the  accuracy  increases  as  the  number  of  Mach  boxes  increase. 
Therefore,  the  maximum  number  of  Mach  boxes  should  be  used.  The  details 
of  the  Transonic  AIC  program  may  be  found  in  Volume  III. 
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3.  The  Supersonic  Box  Method 
n)  The  Technical  Approach  - 

The  basic  work  on  the  supersonic  box  method  was  done  by  Pines, 
Dugundji,  and  Neuringer  (Reference  10)  and  by  Zartarian  and  Hsu  (Reference 
11).  Further  refinements  were  made  by  Moore  and  Andrew  (Reference  12). 
References  7  and  8  also  include  extensions  of  Reference  12  to  interfering 
wing-tail  configurations.  The  fundamentals  of  the  potential  equation  have 
been  presented  above  in  the  introductory  discussion  of  the  technical  approach 
to  the  subsonic  problem  in  Section  A,  and  the  details  of  the  AIC's  and  the 
numerical  problems  of  the  supersonic  box  method  are  discussed  in  the  Unsteady 
Aerodynamics  Report,  Part  III  and  Part  VI,  Section  A,  respectively.  We 
merely  outline  the  basic  procedure  in  this  introductory  summary. 

Pines,  et  al,  (Reference  10)  published  the  first  source  super¬ 
position  method  to  approximate  successfully  the  asrodynmaic  forces  on  an 
oscillating  thin  planar  surface  in  supersonic  flow.  They  employed  Evvard’s 
diaphragm  concept  (Reference  13)  to  handle  subsonic  leading  edges  and  over¬ 
laid  the  surface  and  diaphragm  with  a  grid  of  square  boxes.  For  purposes 
of  calculating  pressures,  they  assumed  that  the  source  strength  over  the 
area  of  each  box  is  a  constant  value  which  satisfies  the  condition  of 
tangential  flow  at  the  center  of  the  box.  The  Mach-oox  procedure  is  basic¬ 
ally  the  same  as  the  method  of  Pines  et  al,  differing  only  in  that  the 
surface  and  diaphragm  is  overlaid  with  a  grid  of  rectangular  boxes,  the 
diagonals  of  which  are  parallel  to  Mach  lines.  As  in  the  Subsonic  and 
Transonic  cases  the  potential  equation  is  written  as  an  integral  equation, 
this  time  relating  the  downwash  to  the  source  strengths,  and  approximated 
by  n  matrix  equation.  The  numerical  difficulties  are  primary  ones  of 
computer  logic  since  the  xones  of  influence  of  a  given  Mach  box  is  limited 
to  the  region  within  the  aft  Mach  lines.  The  matrix  formulation  leads  to 
a  partitioned  form  since  there  are  two  boundary  conditions  to  be  matched. 

The  first  boundary  condition  is  the  downwash  on  the  surface,  and  the  second 
Is  that  there  bo  no  pressure  difference  off  of  the  surface  in  the  diaphragm 
regions.  The  ?,ero  pressure  conditions  leads  to  n  relationship  between  the 
diaphragm  potentials  and  the  surface  potentials,  and  the  surface  downwash 
condition  then  lends  to  the  surface  potentials.  The  surface  pressures  then 
follow  from  the  surface  potentials. 


b)  Program  Capabilities  - 

The  Supersonic  AIC  program  calculates  AIC's  for  systems  with  up  to 
AO  collocation  stations.  The  systems  analyzed  may  be  a  single  surface  or 
tandem  surfaces.  Up  to  20  collocation  stations  per  surface  may  be  used. 

The  following  single  surface  planforms  may  be  analyzed:  rectangular, 
trapezoidal,  delta,  cropped  delta,  and  cropped  trapezoid.  It  is  not 
necessary  to  enter  a  second  surface  into  the  program  when  analyzing  a  single 
surface.  The  trailing  edge  of  the  surface  must  be  straight  and  normal  to 
the  direction  of  flow.  Wien  analyzing  tandem  surfaces,  the  forward  surface 
may  be  any  one  of  the  single  surfaces  previously  defined;  the  aft  surface 
must  be  a  rectangular  surface  with  the  spanwise  dimension  equal  to  the 
spanwise  dimension  of  the  wing. 


The  maximum  allowable  number  of  Mach  boxes  is  U5.  The  layout 
of  tne  boxes  on  the  surface  is  determined  in  the  program,  and  is  a  function 
of  the  number  of  chordwise  boxes  and  the  Mach  number.  The  Mach  boxes  are 
rectangular  with  a  length  in  the  chordwise  direction  equal  to  approximately 
the  root  chord  divided  by  the  number  of  chordwise  Mach  boxes,  and  a  length 
in  the  spanwise  direction  equal  to  the  length  in  the  chordwise  direction 
divided  by  the  square  root  of  the  Mach  number  squared  minus  one.  The  Mach 
boxes  can  be  located  on  the  surface  by  placing  the  center  of  one  box  at 
the  intersection  of  the  root  chord  and  the  trailing  edge  of  the  wing.  It 
is  not  necessary  to  layout  the  Mach  boxes,  as  this  is  done  by  the  program, 
but  this  knowledge  is  useful  in  choosing  the  correct  number  of  chordwise 
M/ich  boxes.  In  general,  the  accuracy  increases  as  the  number  of  Mach 
boxen  increase.  Therefore,  the  maximum  number  of  Mach  boxes  should  be  used. 
The  details  of  the  Supersonic  AIC  program  may  be  found  in  Volume  III. 


PART  IV  -  SAMPLE  PROBLEM 


As  an  example  problem,  a  typical  missile  with  tandem  aerodynamic 
surfaces  is  used.  The  tandem  surfaces  are  in  close  proximity  such  that 
aerodynamic  Interference  between  surfaces  occurs.  Thus,  the  structure  must 
be  modeled  as  one  entity.  The  analysis  will  be  performed  for  the  missile 
free  in  space  with  two  rigid  body  modes  (pitching  and  plunging) ;  thus,  the 
free  boundary  condition  of  the  COFA  program  will  be  employed. 

The  entire  missile  is  flexible.  For  a  missile  free  in  space,  two 
basic  modes  of  flutter  can  exist.  They  are  symmetrical  and  anti- symmetrical 
flutter;  these  modes  of  flutter  are  demonstrated  in  Figure  4.1.  An  analysis 
for  symmetrical  flutter  is  presented.  Due  to  symmetry,  only  one  half  of  the 
structure  needs  to  be  included  in  the  analysis.  A  sketch  of  the  missile  to 
be  analyzed  is  shown  in  Figure  4.2. 

The  first  step  in  the  analysis  is  to  develop  a  structural  idealization 
(mathematical  model)  that  is  consistent  with  the  program  FLUENC.  The 
mathematical  model  established  for  this  example  is  shown  in  Figure  4.3.  The 
structure  was  arbitrarily  fixed  in  translation  at  mass  points  5  and  7,  as  a 
flexibility  matrix  for  a  structure  free  in  space  can  not  be  directly 
calculated.  These  points  will  be  released  later  using  the  free- free  boundary 
condition  of  the  COFA  program.  The  next  step  is  to  run  the  program  FLUENC 
to  generate  the  flexibility  matrix.  The  computer  program  FLUENC  normally 
yields  a  reduced  flexibility  matrix:  i.e.  the  rows  and  columns  pertaining 
to  the  fixed  points  are  eliminated.  This  is  satisfactory  when  the  constrained 
points  are  included  in  a  rigid  component  or  a  fixed  reference  point; 
however,  when  the  constrained  points  are  part  of  a  flexible  structure  or 
component  this  is  not  satisfactory.  A  full  flexibility  matrix  is  required 
In  order  that  the  entire  inertia  of  the  structure  may  be  included.  As  a 
consequence,  the  option  to  punch  the  expanded  flexibility  matrix  for  the 
structure  was  necessary.  This  option  inserts  rows  and  columns  of  zeros  in 
the  flexibility  matrix  corresponding  to  the  constrained  mass  points.  Next 
a  set  of  flight  parameters  were  chosen  and  the  aerodynamics  program  was  run. 
For  tills  example  the  subsonic  AIC  program  was  used,  the  flight  parameters 
used  are 

Much  number  =  .5 

Reduced  Frequency  =  L.Q,  1.5,  2.0 

Altitude  =  sea  level 


If. 


The  computer  output  for  the  subsonic  AIC  program  is  shown  on  pages  38 
through  42  .  As  the  AIC  programs  only  establish  AIC's  for  the  aerodynamic 
surfaces,  the  output  from  the  computer  program  must  be  banded  by  the  proper 
number  of  zeros,  corresponding  to  the  number  of  mass  points  on  the  fuselage. 

The  flexibility  matrix,  the  modified  AIC  matrix,  the  proper  rigid  body  modal 
matrices  corresponding  to  symmetrical  flutter,  and  a  rigid  body  mass  matrix 
corresponding  to  the  arbitrarily  fixed  points  in  the  FLUENC  analysis,  are 
entered  into  the  COFA  program.  The  output  from  the  COFA  program  is  shown 
on  pages  61  through  65  .  The  pertinent  flutter  information  is  shown  on 
page  65  . 

The  data  shown  is  for  the  specific  flight  parameters  entered.  For  the 
output  data  to  be  valid,  the  flutter  velocity  must  equal  the  assumed  input 
velocity.  'Then  this  occurs,  the  damping  value  yielded  is  the  damping 
required  at  that  velocity  and  is  a  point  on  the  flutter  boundary.  Thus, 
determining  a  point  on  the  stability  boundary  is  an  iterative  procedure 
based  upon  varying  the  reduced  frequency  for  a  given  set  of  flight  parameters 
until  the  flutter  velocity  is  approximately  equal  to  the  input  velocity.  In 
general,  this  occurs  at  a  different  reduced  frequency  for  each  flutter  mode; 
thus,  for  one  set  of  flight  parameters,  there  are  as  many  iterative  analyses 
required  as  there  are  modes  requested  in  the  analysis.  The  above  analysis 
determines  one  point  on  the  flutter  boundary.  To  determine  the  entire 
flutter  boundary,  the  above  multiple  iterative  procedure  must  be  repeated 
for  several  sets  of  flight  parameters.  A  typical  "V-g"  curve  is  shown  in 
Figure  4.4  for  a  given  altitude.  Using  a  set  of  these  curves,  for  several 
altitudes,  a  flutter  boundary  of  the  type  specified  in  MIL-A-8870  or 
MIL-M-8856  may  bo  constructed. 

WIhmi  solving  an  iintisymmctricul  flutter  problem  the  technique  is  identical 
to  i  Inn  of  the  nyinmct rical  problem  with  the  exception  of  the  fuselage 
at rucl  in  a|  data  Lnpnt  tor  calculating  the  flexibility  matrix.  For  the 
anil syiraiict  rical  case,  the  fuselage  torsional  stiffness  and  polar  moments 
of  Inertia  rather  than  the  bending  stiffnesses  and  translational  masses 
are  required.  In  addition,  the  fuselage  mass  points  must  be  restrained 
in  translation,  allowing  only  rotation  of  the  fuselage.  The  mass  and 
stiffness  data  for  the  aerodynamic  surfaces  is  identical  to  that  for  the 
syntmel  rical  analysis. 
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As  noted  above,  the  analyses  are  an  iterative  procedure.  Thus,  whether 
an  analysis  is  valid  or  not  is  an  arbitrary  decision,  as  the  output  velocity 
will  never  exactly  be  the  input  velocity.  Normally  when  the  output  velocity 
is  within  five  percent  of  the  input  velocity,  the  results  are  considered 
valid.  For  the  example  problem,  an  input  velocity  of  320  kts  (Mach  .5)  was 
used.  The  output  velocity  for  the  second  mode,  using  a  l/kr  =  .660,  was 
308  kts  (page  65) .  This  normally  would  be  considered  a  valid  point  on  the 
flutter  boundary  for  the  second  mode.  Scanning  the  results  on  page  65  it 
can  be  seen  that  the  valid  data  point  for  the  first  mode  requires  a  l/kf 
value  between  1.00  and  .666;  also  for  the  third  mode  a  value  between  .666 
and  .500.  For  modes  A,  5  and  6,  the  1/k^  value  must  be  greater  than  2.0. 

It  is  noted  that  the  function  of  flutter  velocity  versus  1/k^  is  not 
generally  linear;  thus,  interpolation  will  not  necessarily  yield  the  correct 
1/k^  to  use  in  subsequent  iterations. 
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